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Abstract 

We review two numerical methods related to the Schramm-Loewner evolution (SLE). 
The first simulates SLE itself. More generally, it finds the curve in the half-plane that 
results from the Loewner equation for a given driving function. The second method can 
be thought of as the inverse problem. Given a simple curve in the half-plane it computes 
the driving function in the Loewner equation. This algorithm can be used to test if a 
given random family of curves in the half-plane is SLE by computing the driving process 
for the curves and testing if it is Brownian motion. More generally, this algorithm can 
be used to compute the driving process for random curves that may not be SLE. Most of 
the material presented here has appeared before. Our goal is to give a pedagogic review, 
illustrate some of the practical issues that arise in these computations and discuss some 
open problems. 
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1 Introduction 



This review is about two types of numerical calculations related to the Schramm-Loewner 
evolution (SLE). The first is to simulate SLE itself. More generally, one can consider simulating 
the random curves you obtain in the plane when a random process is used for the driving 
function in the Loewner equation. The second type of simulation is to take a family of random 
curves in the plane and compute the random driving process that generates them through the 
Loewner equation. This is related to SLE since one can test if a given family of random curves 
is SLE by testing if the random driving process is Brownian motion. More generally, it is of 
interest to study the random driving process for random curves that may not be SLE. This 
review is meant to be pedagogic. Most of this material has appeared elsewhere. Our goal is 
to provide the reader with a "how-to" guide that will enable him or her to do state of the art 
simulations related to SLE. 

In the next section we give a heuristic and somewhat atypical introduction to SLE that 
does not involve the Loewner equation. This is followed in section [3] with a quick review of the 
Loewner equation and the usual definition of SLE. The "discretization" of SLE that is used 
in section [2] and discussed in detail in section [3] was studied extensively in p] for a particular 
approximation of the driving function (vertical slits). Reviews of SLE from the mathematics 
point of view include [151 123] and from the physics point of view include [21 El [H] . 

In section H] we consider the numerical algorithm for finding a curve for a given driving 
function in the Loewner equation. Doing this with samples of Brownian motion for the driving 
function gives a simulation of SLE. 

In section [5] we consider the numerical algorithm for finding the driving function for a given 
curve. One motivation for doing this is that it gives a way to test if a given model is SLE by 
testing if the driving process is Brownian motion. Several works have considered models for 
which the connection with SLE is not clear, including domain walls in spin glasses [H [7] and 
turbulence [5l [6]. Another motivation is to study the driving process for massive scaling limits 
of off-critical models [1 [3 [H] 

Both of the numerical algorithms we study are closely related to the zipper algorithm |T^[TB]. 
This is an algorithm for numerically finding the conformal map of a given simply connected 
domain onto a standard domain such as the unit disc. Much of the work described in this 
review grew out of conversations with Don Marshall and Stephen Rohde. 

2 An introduction to SLE 

In this section we will give a heuristic introduction to SLE. The standard definition of SLE 
uses the Loewner equation from complex analysis. We will give a different definition of the 
process that does not use the Loewner equation. This view of SLE is well known, but is not 
typically discussed in reviews of SLE. The approach to SLE that we present is closely related 
to the numerical algorithms we will discuss. In the next section we will see how this approach 
is related to the usual definition using the Loewner equation. 
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Let EI denote the upper half of the complex plane, 

m= {z : Im{z) > 0} (1) 

Fix an angle 6 G (0,7r/2] and a length p > 0. Let f+{z) be the conformal map which takes H 
onto EI \ {re*^ : < r < p}, the upper half plane minus the line segment from to pe*^. This 
map is not unique. We make the choice unique by requiring 

/+(oo) = oo 

/;(oo) = 1 

/+(0) = pe'' (2) 
The first two conditions mean that the Laurent series of /+ about oo is of the form 

f^(z) = z + Co + - + ^ + --- 
z z-^ 

(For the reader familiar with the Loewner equation, we note that this is not the "hydrodynamic" 
normalization which would require that Cq = in the Laurent expansion instead of the third 
condition in ([2]).) The map /+ is illustrated by the upper left picture in figured! The grid 
shown is the image under the conformal map of the uniform rectangular grid in the upper half 
plane. Let /-(z) be the analogous conformal map for the segment from to pe*^'^"^-'. (So the 
range of /_ is the reflection of the range of /+ about the vertical axis.) 

Consider composing two of these maps, e.g., /+ o /_. The effect of the second map in the 
composition will be to push the line segment created by the first map into the upper half plane 
and bend it somewhat. Because we have required that these maps send to the tip of the line 
segment, the lower endpoint of the image of the first slit under the second map will be the tip 
of the second slit. In other words the image of EI under the composition will be HI with a curve 
removed. The map /+ o /_ is illustrated by the picture in the upper right of figure [H 

We can compose multiple copies of /_ and /+ and the resulting conformal map will send 
the half plane onto the half plane minus a curve. We choose the maps randomly. Let 
be a sequence of independent, identically distributed random variables with = ±1 with 
probability 1/2. For positive integers n consider the conformal map 

Fn = fxi o fx2° fx3° ■ ■ ■ ° /x„ (3) 

(There is a slight abuse of notation here: f±i means f±.) The picture in the lower left of figure 
[1] illustrates an example of F^, and the picture in the lower right an example of Fiq. 

The conformal map F„ will map HI onto EI \ 7„ where 7„ is a curve in the upper half plane 
starting at 0. Because of the order of the X'^s in ([3]), the curve 7„+i will be an extension of the 
curve 7„. So we can let oo to get an infinite curve 7. The SLE curve is then obtained by 
taking the scaling limit p ^ 0. The angle a is related to the usual parameter k, for SLE. Using 
results that will appear latter, one can show that the relation is 

4(2a-l)2 

K = — 

a{l — a) 
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Figure 1: The figures illustrate the random composition of a sequence from the maps 
/_ and /+. The numbers of maps in the compositions are 1, 2, 5 and 10. 

When we use a piecewise smooth approximation to the driving function in the Loewner 
equation, the curve 7 will be simple (non-intersecting). It is a subtle question whether the 
curves one obtains in the limit p — > are simple. For k < A, SLE produces a simple curve [20], 
and it is natural to conjecture that 7 converges to this SLE curve. We do not prove this, and 
we are not aware of any proof in the literature. For k > A the random set produced by SLE 
is not even a curve |20]. It is generated by a non-simple curve, called the SLE trace, in the 
sense that the SLE set at time t is the complement of the unbounded connected component 
of the half plane minus the curve up to time t. It is natural to conjecture that 7 converges in 
distribution to the SLE trace, but again we do not prove this and are not aware of any proof 
in the literature. Closely related questions are considered in [2]. 
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3 The Loewner equation 



We will now quickly review the Loewner equation from complex analysis and see how it is 
related to the definition of SLE that we gave in the previous section. The Loewner equation 
provides a means for encoding curves in the upper half plane that do not intersect themselves 
by a real-valued function. In fact, it applies to more general growth processes in the half plane, 
but for the moment we restrict our attention to curves. Let ^{t) be a simple curve which lies 
in H for < t < cxD and starts at the origin, i.e., 7(0) = 0. Let 7[0, t] denote the image of 7 up 
to time t. Then HI \ 7[0, t] is a simply connected domain. So there is a conformal map gt from 
this domain to H. This map is not unique. We choose the map that satisfies 

g,(z)=Z + ^ + 0{^), Z^OO (4) 

z \z\^ 

The coefficient C{t) is called the half-plane capacity of 7[0,t]. It is known to be increasing in 
t, so we can parametrize the curve so that C{t) = 2t. Then gt satisfies Loewner's differential 
equation 

dgtiz) 2 

= —r\ F7' 9o[z) = z (5) 

dt gt{z) - Ut 

for some real valued function Ut on [0, cxd). This statement is not obvious, and we refer the 
reader to [15] for a proof. The function Ut is often called the driving function. We emphasize 
that while gt{z) is complex valued, the driving function Ut is real- valued. 

Note that gt goes in the opposite direction of the maps in the previous section, i.e., it sends 
the half plane with a curve deleted onto the half plane while the previous maps sent the half 
plane onto the half plane minus a curve. We should also note that gt is normalized differently 
since the constant term in (jlj) vanishes. So gt{l{t)) is not the origin. In fact it is Ut- (To be 
precise, gti^jit)) is not defined since •yit) is on the boundary of the domain of gt- Its image 
under gt must be defined by a limiting process.) 

If our simple curve in the half plane is random, then the driving function Ut is a stochastic 
process. Schramm's wonderful discovery was that if the scaling limit of a two-dimensional model 
is conformally invariant and satisfies a property usually called the domain Markov property, 
then this stochastic driving process must be a Brownian motion with mean zero [21]. The only 
thing that is not determined is the variance. Schramm named this process stochastic Loewner 
evolution or SLE; it is now often referred to as Schramm-Loewner evolution. 

The solution to ([5]) need not exist for all times t since the denominator can go to zero. We 
let Kt be the set of points z in IH for which the solution to this equation no longer exists at 
time t. If we start with a simple curve and define gt as we did above, then Kt will be 7[0,t]. 
But if we start with a continuous driving function Ut and solve the Loewner equation, Kt will 
only be a curve for sufficiently nice Ut. (Just what sufficiently nice means is a subtle question 
[T7].) For other Ut, Kt can be a more complicated growing set. In particular, when Ut is a 
Brownian motion, Kt may not be a curve. In our simulations, even in the cases where Ut is not 
sufficiently nice, our approximation to Ut will be nice enough that it produces a curve. So in 
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the following we will always take Kt to be a curve, but the reader should keep in mind that in 
some cases this curve is approximating a more complicated set. 

Let t,s > 0. The map gt+s niaps E[\7[0,t + s] onto H. We can do this in two steps. We first 
apply the map gs- This maps H[\7[0, s] onto H, and it maps EI\7[0, t + s] onto EI\5fs(7[s, t + s]). 
Let cjt be the conformal map that maps HI \ gs{'~f[s, t + s]) onto EI with the usual hydrodynamic 
normalization. Then cjt o gs will map EI \ 7[0, t + s] onto EI and satisfy (jl]). There is only one 
such conformal map, so 

9s+t = gt0 9s, i-e., gt = gs+t0 9s^ (6) 

It we think of s as being fixed and t as the time variable, then the function gt is also a solution 
of the Loewner equation 

^9t{z) = ^Qs+t o 9^^{z) = 3^^— -— = _ , , ^ (7) 

at dt gs+t og^\z) ~ Us+t 9t[z) - Us+t 

and satisfies go{z) = z. Thus gt{z) is obtained by solving the Loewner equation with driving 
function Ut = Us+t- This driving function starts at Us, and so the curve associated with gt 
starts at Us- 

We now introduce a partition of the time interval [0, oo): = to < < ^2 < ■ ■ ■ < • ■ ■, 
and define 

9k = 9t, o 9^,\ (8) 

So 

9tk = 9k° 9k-i o gk-2 ° ■ ■ ■ °g2°gi (9) 

By the remarks above, gk is obtained by solving the Loewner equation with driving function 
Ut^_^+t for t = to t = Afe, where = — tk-i- The image of EI under gk is EI minus a "cut" 
starting at Ut^^^- So if we shift it by defining 

gk{z) = -gk{z + Ut,_,)-Ut,_,, (10) 

then gk is obtained by solving the Loewner equation with driving function Ut^_^+t — Ut^_^ for 
t = to t = Afc. This driving function starts at and ends at 5k where 5k = Ut,^ — Ut^._^. So 
this conformal map takes EI minus a cut starting at the origin onto H. The inverse of this map, 

9k\^) = 9k\z + Ut,_,) - Ut,_„ (11) 

takes EI and introduces a cut which begins at the origin. 

There are two general types of simulations we would like to do. Given a driving function we 
want to find the curve it generates. And given a curve we want to find the corresponding driving 
function. For both problems the key idea is the same. We approximate the driving function on 
the interval [tk-i,tk] by a function for which the Loewner equation may be explicitly solved. 
So the maps gk and gk can be found explicitly. Eq. ([9]) can then be used to approximate gt. 
We will consider two explicit solutions of the Loewner equation which we will refer to as "tilted 
slits" and "vertical slits." 
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For tilted slits, let xi,Xr > and < a < 1. Then define 



Then / maps EI to HI \ F where F is a line segment from to a point pe*"'^. The length p can be 
expressed in terms of xi,Xr and a. This map sends [—xi,Xr] onto F. Unfortunately, its inverse 
cannot be explicitly computed. For the inverse to satisfy the normalization (jlj), we must have 

(1 — a)xi = axr (12) 

Straightforward calculation shows if we let 
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1 — a J \ \ a 

then it produces a slit with capacity 2t. We know that gt = f^^ must satisfy the Loewner 
equation ([5]) for some driving function Ut- More calculation shows that the driving function is 

r- 1 — 2a 

Ut = CaVt, Ca = 2 (13) 

y «(1 — a) 

The change in the driving function over the time interval [0, A] is 

5 = c„v^ (14) 

The original map had three real degrees of freedom, a,xi,Xr- The condition ( |T2l) reduces this 
to two real degrees of freedom, a and t. So if we are given S and A or given p and a, then the 
map is completely determined. 

Vertical slits correspond to an even simpler solution of the Loewner equation. Let 
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(z) = ^/iz-sy + u + s 



Then it is easy to check that gt satisfies Loewner's equation with a constant driving function, 
Ut = S. Since the driving function does not start at 0, the curve will not start at the origin. The 
curve is just a vertical slit from 6 to 6 + 2i\ft. Using vertical slits means that we approximate 
the driving function by a discontinuous piecewise constant function. This will produce a Kt 
which is not a curve. 

Our numerical studies only use tilted slits and vertical slits for the explicit solutions for the 
Loewner equation. Another possibility is to use a linear driving function. If we let ht = gt — Ut, 
then the differential equation for ht can be solved by separation of variables. The solution is 
not completely explicit - it must be expressed in terms of a function that is defined implicitly 
by a transcendental equation. 
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4 From the driving function to the curve 



Our primary motivation is to simulate SLE, i.e., to compute the curve when the driving function 
is Brownian motion. But our discussion is more general, and the following algorithm can be 
used to calculate the curve corresponding to any driving function Ut- 

There are a variety of conformal maps that occur in this paper, and we have denoted them 
by letters that indicate what they do. Maps denoted with g are solutions of the Loewner 
equation with a driving function that starts at 0. So they map the half plane minus a curve 
starting at the origin onto the half plane, sending the tip of the curve to the final value of the 
driving function. We use g for solutions to the Loewner equation when the driving function 
does not start at 0. In this case the curve starts at the initial value of the driving function and 
the map still sends the tip to the final value of the driving function. If we follow a map g hj a. 
real translation that takes the final value of the driving function to 0, we get a map that takes 
the half plane minus a curve onto the half plane and sends the tip to the origin. We denote 
such maps by h. (Note that such maps do not satisfy the Loewner equation.) Finally, we use 
/ to denote maps that are inverses of maps h. So they take the half plane onto the half plane 
minus a curve and sends the origin to the tip. 

Let = to < ^1 < ^2 < ■ ■ ■ < i^n be a partition of the time interval [0,t]. The SLE curve 
is given by •yit) = g^^iUt). Let = g^^iUt^.). We will only consider the points on this 
curve which correspond to times t = t^. One could consider other points on the curve, but the 
distance between consecutive Zk is already of the order of the error in our approximation, so 
there is no reason to consider more points. By ([9]) the points z^ are given by 

Zk = gi^ og2^° gk-i ° 9k\Ut,) (15) 

Recall that if we solve the Loewner equation with driving function Ut,._^+t — Ut^_^ for t = to 
t = A/;, the result is gk{z) where 

gk{z) = gk{z + Ut^_^) - Ut^_^ (16) 

Define 

hk{z) = gkiz) -5k = g{z + Ut,_J - Ut, (17) 
where Sk = Ut, — Ut,_^. Then 

hk o hk-i o o hi{zk) =gk° gk~i o o gi{zk) - Ut, = 0. (18) 

Let 

fk = (19) 

So 

^fe = /lO/2O---o/,(0) (20) 

As noted before, gk maps EI minus a small curve onto H. The driving function ends at 6k, 
so gk sends the tip of the curve to 6k- It follows that hk{z) = gk{z) — 6k maps H minus the 
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small curve onto H and sends the tip to the origin. So fk = maps EI onto EI minus the 
small curve and sends the origin to the tip of the curve. Thus the functions fk are analogous 
to the functions f± from section [2] in that they all introduce a small cut into the upper plane 
and send the origin to the tip of the cut. Note the similarity of ( l20l) to ([3]). 

As discussed before, we define Ut on each time interval tk-i < t < tk so that gk{z) may be 
explicitly computed. There are two constraints on g^. The curve must have capacity 2Ak and 
Qk must map the tip of the curve to 6k- Any simple curve satisfying these two constraints and 
starting at the origin will correspond to a solution of the Loewner equation for some driving 
function which goes from to 6k over the time interval [0, Afc]. So our approximation can be 
thought of as replacing the driving function by a new driving function that agrees with the 
original one at the times tk but differs in between those times. 

Different choices of how we define Ut on each time interval give us different discretizations. 
As we will see, this choice will not have a significant effect. Of much greater importance is how 
we choose the A^, and 6k- 



Figure 2: SLE with k = 8/3 with fixed At. There are 10, 000 points. 

If we want to simulate SLE, the 6k should be chosen so that the stochastic process Ut 
will converge to ^/k times Brownian motion as — * oo. One choice is take the 6k to be 
independent normal random variables with mean zero and variance nAk- If we do this, then Ut 
and y/nBt will have the same distributions if we only consider the times tk- Another possibility 
is to approximate the Brownian motion by a simple random walk. This is done by using a 
uniform partition of the time interval and taking the 6k to be independent random variables 
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with 5k = iy^/tAfc where the choices of + and — both have probabihty 1/2. This is what we 
were doing in section [2l 

The simplest choice for is to use a uniform partition of the time intervaL For values of k 
which are not too large this works reasonably well. Figure [2] shows a simulation using /t = 8/3 
with 10, 000 equally spaced time intervals. However, for larger values of k, uniform are a 
disaster. Figure [3] shows a simulation with k = 6 and 10, 000 equally spaced time intervals. 
Clearly something has gone wrong. To see just how badly wrong things have gone the reader 
should compare this figure with figure H] which uses the same sample of Brownian motion. 




Figure 3: SLE with k = 6 with fixed At. There are 10, 000 points. 

To understand the effect seen in figure [3] we give an equivalent definition of the half plane 
capacity C of a set A. We originally defined it by 

g{z) = z + - + 0{^) 

where g maps EI \ A onto H. A more intuitive definition is 

C = lim yE'y[Im{Br)] 

where Bt is two-dimensional Brownian motion started at iy. The stopping time r is the first 
time the Brownian motion hits A or R. From the point of view of this two-dimensional Brownian 
motion, parts of the curve can be well hidden by earlier parts of the curve and so have very 
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Figure 4: SLE with k = 6 with adaptive At. There are 35, 000 points. 



httle capacity. So what looks hke a "long" section of the curve has very little capacity and so 
gets very few points approximating it. 

To do better we will use non- uniform A^. In fact they will depend on the sample of the 
Brownian motion and so we refer to this method as "adaptive A^." (I learned this idea from 
Stephen Rohde [19j.) Fix a spatial scale e > 0. We start with a uniform partition of the time 
and compute the points along the curve. Then we look for points Zk such that \zk — Zk-i\ > e. 
For these time intervals [tk-i,tk], divide the interval into two equal intervals. We then sample 
the Brownian motion at the midpoint of [tk-i,tk] using a Brownian bridge. (This just means 
that to choose the value of the driving function at the midpoint of [tk-i, tk] we use a Brownian 
motion conditioned on the values we already have for it at t^-i and t^.) Then we recompute all 
the Zk. (There will of course be more of them than before.) Note that we must recompute all 
the points since even at times which appeared in the time partition before, the corresponding 
point on the curve will change. We repeat this until we have \zk — Zk-i\ < e for all k. 

Our approximation can be thought of as approximating the driving function by a concate- 
nation of driving functions on short time intervals for which the Loewner equation is exactly 
solvable. It is important to consider the effect of the choice of which exactly solvable driving 
functions we use. To do this we compare the curves we get using tilted slits for the elementary 
maps with the curves we get using vertical slits. We carry out the adaptive simulation just 
described using tilted slits. We then use the same A^ and 6k, i.e., the same partition of the time 
interval and the same sample of Brownian motion, but with vertical slits. For n = 8/3, figure 
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[5] shows the tilted shts curve vs. the vertical slits curve. The vertical slits do not produce a 
curve. What we have plotted is the following. We compute the points Zk and then just connect 
them with a straight line. In figure [5] it is almost impossible to distinguish the two curves. An 
enlargement of part of the curves is shown in the inset. Even in the enlargement the difference 
is quite small. The curves have a relatively small number of points (about 6,000), and in the 
enlargement we have plotted the points for the tilted slit curve. The difference between the 
two curves is on the order of the distance between these points. 

Figure [6] shows the same thing with k = Q. In the enlargement one can see deviations 
between the two curves, but the size of the deviations is again on the same scale as the distance 
between adjacent points on the curve. 

It is interesting to note that there is what one might call a stability to the approximation 
we are using. The difference between the two curves in figures |5] and [6] fiuctuates with time, 
but it does not grow with time. In other words, the errors from approximating the true driving 
function over the short time intervals do not appear to accumulate. 



tilted slits 








Figure 5: A comparison of the curves obtained using tilted slit maps and vertical slit 
maps with k = 8/3. 
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Figure 6: A comparison of the curves obtained using tilted slit maps and vertical slit 
maps with k = 6. 
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5 From the curve to the driving function 



We now consider what one might call the inverse problem. Given a simple curve 7, we want to 
compute the corresponding driving function. 

Let 7(5) be a parametrized simple curve in H. In almost all applications, the parametrization 
of the curve is not the parametrization by capacity. Let Qs be the conformal map which takes 
H \ 7(0, s] onto H, normalized so that for large z 

gs{z) = z+^ + 0{\), (21) 
z z 

The coefficient C{s) is the half-plane capacity of 7[0,s]. The value of the driving function 
t — C{s) /2 is Ut — gs{j{s)). Thus computing the driving function essentially reduces to 
computing this uniformizing conformal map. 

Let Zo, Zi, ■ ■ ■ , Zn be points along the curve 7 with zq = 0. In many applications these 
are lattice sites. We will find a sequence of conformal maps hi, i = 1,2, such that 

hk o hk-i o ■ ■ • o hi{zk) ~ 0. Then o h^^i o ■ ■ ■ o hi sends H \ 7 to H where 7 is some curve 
that passes through Zq,Zi, - • ■ z^ and so approximates 7. Suppose that the conformal maps 
hi,h2, • ■ • ,hk have been defined with these properties. Let 

Wk+i = hk o hk-i o ■ ■ ■ o hi{zk+i) (22) 

Then w^+i is close to the origin. We define h^+i to be a conformal map that sends EI \ 'jk+i to 
EI where 7,^+1 is a short simple curve from to Wk+i- We also require that /i^+i sends w^+i to 
the origin. As before we choose the curve 7^+1 so that hk+i is explicitly known; possible choices 
include "tilted slits" and "vertical slits." Note that for both of these maps there were two real 
degrees of freedom. They will be determined by the condition that hk+i{wk+i) = 0. 

Let 2Aj be the capacity of the map hi, and Si the final value of the driving function for hi. 

So 

2A- 1 

hi{z) = z-6i + —^ + 0{-) (23) 



z^ 



Then 



where 



2t 1 

hk o hk-i o • • • o hi{z) ^z-Ut + — + 0{—) (24) 



z ' z"^' 



i=l 

k 



Ut^Y^^i (26) 

i=l 

Thus the driving function of the curve is obtained by concatenating the driving functions of 
the elementary conformal maps h^. 
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6 Faster algorithms 



In this section we show how to speed up both the algorithm for computing the curve 7 given 
the driving function Ut and the algorithm for computing the driving function Ut given a curve 
7. We start with the first algorithm. One of the main motivations is a fast algorithm for 
simulating SLE, but our fast algorithm is applicable to other driving functions as well. 

Recall that points on the approximation to the SLE trace or more generally the curve 7 are 
given by eq. ( |20l) which says 

Zk = fiof20---o /,(0) (27) 

The number of operations needed to compute a single Zk is proportional to k. So to compute 
all the points Zk with k = 1,2, - ■■ N requires a time 0{N'^). The computation of Zk does not 
depend on any of the other Zj. Depending on what we want to compute, we may only need to 
compute a subset of the points Zk- (For example, if we are only interested in zj^ = 7(tAr), the 
time required is 0{N) not 0{N^).) For a typical point Zk, the time to compute it is 0{N) for 
the above algorithm. Our goal is to develop an algorithm for which this time is 0{N^) with 
p <l. 

Our algorithm groups the functions in fl27j) into blocks. We denote the number of functions 
in a block by h. Let 

= ° f(j-i)b+2 ° fjb (28) 

If we write k as k = mb + / with < / < 6, then we have 

2;fe = Fi O F2 O ■ ■ ■ O F„ O /^b+i O o ■ ■ ■ o fmb+liO) (29) 

The number of compositions in ( l29|) is smaller than the number in ( !27l) by roughly a factor of 
6 if 6 is smaller than m, i.e., if k is bigger than 6^. 

Unfortunately, even though the fi are explicit and relatively simple, the Fj cannot be 
explicitly computed. Our strategy is to approximate the fi by functions whose compositions 
can be explicitly computed to give an explicit approximation to Fj. For large z, fi{z) is given 
by its Laurent series about oo. One could approximate fi by truncating this Laurent series. 
Our approximation is of this nature, but slightly different. 

Let 7 : [0, t] — i> H be a simple curve in the upper half plane with 7(0) = 0. Let f{z) be the 
conformal map from HI onto EI\7[0, t]. We assume that / is normalized is the same way as our 
fi, i.e., f{oo) = 00, f'{oo) = 1 and /(O) = 7(t). Let a, 6 > be such that [— a, &] is mapped 
onto the slit 7[0,t]. Then / is real valued on {—00,— a] U [b,oo), and so / has an analytic 
continuation to C \ [—a, b] by the Schwartz reflection principle. We denote this extension by 
just /. 

Let R = max{a,6}, so / is analytic on {z : \z\ > R} and maps 00 to itself. Thus f{l/z) is 
analytic on {z : < \z\ < 1/-R}. Since our assumptions on / imply it has a simple pole at the 
origin with residue 1, we have 

00 

f{l/z) = l/z + J2ckz' (30) 

k=0 
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This gives the Laurent series of / about oo. 



f{z) = z + J2ckZ-' (31) 



k=0 



This Laurent series is a natural approximation to use for / when z is large. However, we will 
use a different but closely related representation. 

Define f{z) = l/f{l/z). Since f{z) does not vanish on {\z\ > R}, f{z) is analytic in 
{z : \z\ < 1/-R}. Our assumptions on / imply that /(O) = and /'(O) = 1. So / has a power 



series 



f{z)=z + J2(^j^' (32) 
i=2 

The radius of convergence of this power series is easily shown to be 1/R. Note that the 
coefficients of the power series of / are the coefficients of the Laurent series of 1//. 

The primary advantage of our power series over the Laurent series is its behavior with 
respect to composition. 

Thus 

(/l O = A O /2 (34) 

Our approximation for f{z) is to approximate f{z) by the truncation of its power series at 
order n. So 

1 



.i=o 



z-^ 



(35) 



For each fi we compute the power series of fi to order n. Using these and (IMIl . we compute 
the power series of Fj to order n. Let l/Rj be the radius of convergence for the power series 
of Fj. Now consider evaluating the composition in equation fl29|) . Let z be the argument to 
Fj. If z is large compared to Rj, then Fj{z) is well approximated using the power series of Fj. 
We introduce a parameter L > 1 and use the power series of Fj to compute Fj{z) whenever 



\z 



> LRj. When \z\ < LRj, we just use ( 128|) to compute Fj{z). The argument of Fj is random, 
and so whether or not we can approximate a particular Fj using these power series is random. 
As part of the algorithm we must compute Rj. This is easy. Rj is the smallest positive number 
such that Fj{Rj) and Fj{—Rj) are both real. 

In addition to the choice of simple curves we use (tilted slits, vertical slits, ....), there are 
three parameters in our algorithm, b is the number of functions composed in a block, n is the 
order at which we truncate our series approximation. L is the scale that determines when we 
use series for Fj. The parameter b has little effect on the accuracy of the algorithm and we 
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should choose it to make the algorithm run as quickly as possible. Eq. (!29|) suggests that b 
should vary with N as and experiments bear this out. 

The choice of n involves a tradeoff of speed vs. accuracy. Larger n means more terms in 
the series, hence slower but more accurate computations. We typically use n = 12. 

The parameter L will determine how fast the series converges. Roughly speaking, the series 
will converges at least as fast as the geometric series Yin • "^^^ choice of L also involves 
a tradeoff of speed vs. accuracy. Larger L means the series converges faster and so is more 
accurate. But it also means that we use the block functions Fj less frequently, and so the 
computation is slower. We typically use L = 4. 

A detailed study of the effects of the choices of b, n and L can be found in [12]. This paper 
also studies the time to compute a point on the curve and finds it is 0{N^) with p approximately 
0.4. To illustrate the accuracy of our series approximation we compute an SLE curve for k = 6 
with and without the series approximation. We use the same Brownian motion sample path for 
both curves. We typically take n = 12 and L = 4. With these choices the difference between 
the curves obtained with and without the series approximation is extremely small and cannot 
be seen in plots of the curves. If we reduce n to only 6 we can begin to see the effect of the series 
approximation. Figure [7] shows the two curves we get for k = 6 and the same sample of the 
driving process when we use n = 6. One can only distinguish the difference in the enlargement 
and even then it is small. 




Figure 7: Two curves for SLE with k = 6 are shown. They use the same Brownian motion 
sample path but one uses the series approximation and the other does not. 
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We now consider the algorithm for computing the driving function of a given curve. The 
number of operations needed to compute a single Wk+i is proportional to k. So to compute all 
the points Wk+i, and hence the approximation to the driving function, requires a time 0{N'^). 
The idea for improving this is the same as before - we group the functions we are composing 
into blocks and approximate the composition F of the functions in a block using the power 
series of F. The only minor difference is that the order of the conformal maps in fl22|) is the 
opposite of that in fl20l) . We continue to denote the number of functions in a block by b. Let 

Hj = hjh o hjf,_i o ■ • ■ o /i(j_i)b+2 o ^(i-i)b+i (36) 
If we write k as k = mb + r with < r < 6, then fl22p becomes 

U!k+l = hmb+r O hmb+r-1 O ■ ■ ■ O hmb+1 ° Hm O O ■ • ■ O Hi{Zk+l) (37) 

As before, the hi are relatively simple, but the composition Hj cannot be explicitly computed. 
We approximate hi by the power series of hi and compute the approximations to the composi- 
tions in (!36l) just once rather than every time we compute a Wk- 

Recall that hi is normalized so that hi{oo) = oo and hi{oo) = 1. It maps H minus a simple 
curve near the origin to H, sending the tip of the curve to the origin. Let h denote such a 
conformal map. Let R be the largest distance from the origin to a point on the curve. Then h 
is analytic on {z G HI : \z\ > R}. Since h is real valued on the real axis, the Schwarz reflection 
principle says it may be analytically continued to {z E C : \z\ > R}. Moreover, it does not 
vanish on this domain. We could approximate h by its Laurent series about oo, but as with 
the first algorithm it is better to use the power series of h{z) = l/h{l/z). Note that the radius 
of convergence of this power series is 1/R. 

As before, the advantage of working with the power series of h is its behavior with respect 
to composition: {hi o /i2)" = hi o h2 Our approximation for hi{z) is to replace hi{z) by the 
truncation of its power series at order n as we did in eq. (l35l) . The approximation of hi and of 
Hj proceeds as for the first algorithm. For each hi we compute the power series of hi to order 
n. We then use them to compute the power series of Hj to order n. As before we introduce 
a parameter L > 0. Let 1/Rj be the radius of convergence for the power series of Hj. Now 
consider equation fl37|) . If the argument z of Hj satisfies 1^1 > LRj, then we approximate Hj[z) 
using the power series of Hj. Otherwise we just use to compute Hj{z). The argument of 
Hj is random, so as before whether or not we can approximate a particular Hj by its series is 
random. 

We need to compute Rj. Consider the images of ZQ^i)b, ■ ■ ■ Zjb-i under the map 

Hj_ioHj_20- ■ -oHi. The domain of the conformal map Hj is the half-plane EI minus some curve 
Tj which passes through the images of these points. The radius Rj is the maximal distance 
from the origin to a point on Tj. This distance should be very close to the maximum distance 
from the origin to images of -2(j-i)f,+i, ■ • • Zjb-i under Hj^i o Hj^2 o ■ ■ • o Hi. So in our 

algorithm we approximate Rj by the maximum of these distances. 

To compute the driving function without using the power series we must compute all the 
points Wk. So if we do not use the power series, the time needed is 0{N'^). The improvement 
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in the speed of the algorithm from using the power series approximation is studied in |13j . 
Numerical experiments indicate it is 0{N^) with p approximately equal to 1.35. 

7 Conclusions and open problems 

We have reviewed numerical methods for taking a driving function and finding the curve pro- 
duced by the Loewner equation and for taking a curve in the half plane and finding the cor- 
responding driving function. Both methods are based on approximating the driving function 
over short time intervals by a function for which the Loewner equation may be solved explicitly. 
The solution of the Loewner equation over the entire interval is then given by a composition 
of such maps. Our numerical studies used as the simple maps the conformal maps that pro- 
duce a vertical slit or a tilted slit in the half plane. The difference in the results when we use 
vertical slits or tilted slits is small. The vertical slit map is considerably faster and simpler to 
implement, so we see no reason to use the tilted slit map. To simulate SLE effectively it is 
imperative that the choice of time intervals be done in a way that depends on the sample of 
the driving function so that sections of the curve that correspond to small changes in capacity 
are computed accurately. 

The speed of both algorithm can be greatly increased by using power series approximations 
of certain analytic functions. This approximation is quite accurate and the errors from it are 
insignificant compared to the effect of changing the number of points used on the curve or 
compared to the difference between using vertical slits or tilted slits in the algorithm. 

We end with a discussion of a variety of open problems related to these two algorithms. 

We have only discussed the simulation of chordal SLE. In chordal SLE the random curve 
goes between two boundary points, e.g., the origin and infinity in the half plane. In radial SLE 
the random curve goes between a boundary point and an interior point, e.g., the point 1 and 
the origin in the unit disc. The simulation of radial SLE is similar. Can one use the ideas we 
used to speed up the simulation of chordal SLE to speed up the simulation of radial SLE? 

Instead of taking the scaling limit at the critical point, one can consider off critical models 
and take the scaling limit in such a way that it has a finite correlation length. What can you 
say about the driving process for this scaling limit ? For percolation it is know to be rather 
nasty [18]. See also [118]. 

There are several methods for numerically computing the conformal map of a given simply 
connected domain onto a standard domain, like the unit disc. One of these methods, the zipper 
algorithm [Ml [T6] , reduces the problem to that of finding the conformal map from the half plane 
minus a curve to the half plane. So the power series approximation that we use also provides 
a faster version of this algorithm. How does this faster version compare to other methods for 
finding the conformal map from a simply connected domain to a standard domain [lOl I22i| ? 

As discussed in section |2l it is natural to conjecture that the discrete SLE curve 7 introduced 
in that section converges to the SLE curve for k < 4 and converges to the SLE trace which 
generates the SLE hull for k > 4. Prove this. Part of the problem is figure out the sense in 
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which they converge. 

For the inverse problem of finding the driving function for a given curve, there is an analogous 
convergence question. Show that as the number of points used to approximate the curve goes 
to infinity, the computed driving function converges to the true driving function. Marshall and 
Rohde have proved convergence for a particular variant of the zipper algorithm p^. 

As discussed in section |U there is a certain stability to our approximation of the curve 
generated by a given driving function. Explain this stability. 
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